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In this lecture we argue that the fluctuations of Dirac eigenvalues on the finest 
scale, i.e. on the scale of the average level spacing do not depend on the underlying 
dynamics and can be obtained from a chiral random matrix theory with the same 
low energy effective theory. We present three pieces of evidence supporting that 
such microscopic correlations of lattice QCD Dirac spectra are given by chiral 
random matrix theory. First, we find that the spectral correlations of eigenvalues 
in the bulk of the spectrum obey the Dyson-Mehta-Wigner statistics. Second, we 
show that the valence quark mass dependence for sufficiently small quark masses, 
as calculated by the Columbia group, can be obtained from the microscopic spectral 
density of chiral random matrix theory. Third, in the framework of chiral random 
matrix models, we present results showing that the microscopic spectral density is 
strongly universal, i.e. is insensitive to the details of the probability distribution 



1 Introduction 

Although not a physical observable itself, the spectrum of the Euclidean QCD 
Dirac operator is an essential ingredient for the calculation of hadronic cor- 
relation functions. It also serves as the order parameter for the chiral phase 
transition. This becomes clear by expressing the quark propagator in terms of 
the eigenfunctions and the eigenvalues of the Dirac operator 



where 

ijD(j) k = A fc (/)fc, (2) 

and i'yD is the Euclidean Dirac operator. Because {75,17!)} = 0, all nonzero 
eigenvalues occur in pairs ±Afe. The eigenvalues near zero are of great impor- 
tance for the propagation of light quarks. This is the reason that instanton 
field configurations, with one zero eigenvalue per instanton, determine the main 
characteristics of hadronic correlation functions El. 
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With regards to chiral symmetry breaking, the Dirac operator is directly 
related to the chiral condensate 
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which is the celebrated Banks-Casher formulaa. Here, the spectral density is 
defined as 



and the average is over gauge field configurations weighted according to the 
Euclidean action. The space-time volume is denoted by V. The limit is defined 
such that the thermodynamic limit is taken before the chiral limit. 

In this lecture, I wish to argue that some properties of the Dirac spectrum 
are completely determined by the global symmetries of the QCD partition 
function and are not sensitive to the dynamics of the theory. In general such 
universal properties are given by fluctuations of the eigenvalues on the finest 
scale, i.e. on the scale of individual level spacings. Because the spectrum is 
symmetric about zero, we have to distinguish spectral correlations near zero 
virtuality from correlations in the bulk of the spectrum. Although the overall 
spectral density is certainly not universal, we wish the convince the reader 
that the sperlcal density near zero, on the scale of the average level spacing, 
is univcrsalU'u'Q. It is called the microscopic spectral density. In the bulk of 
the spectrum, the spectral correlations will be shown to obey universal level 
statisticsu. From the study s£ quantum chaos, we know that the latter correla- 
tions are strongly universal Era. In the framework of a random matrix model we 
will show thatr-the microscopic spectral density has strong universality prop- 
erties as wellBll3. Finally, we make contact with lattice QCD calculations. In 
this context the Dirac spectrurn_svas investigated via the valence quark mass 
dependence of chiral condensates, and, for a relatively small number of config- 
urations, by a full diagonalization of the Dirac operator Both the valence 
quark mass dependence LL-3 and the eigenvalue correlations in the bulk of the 
spectrumQ are shown to be in complete agreement with random matrix theory 
(RMT) within its range of applicability. 




(4) 



k 



2 



2 Universality 

Leutwyler and SmilgaEl have argued that in the range 
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where L is the linear size of the Euclidean box, m the quark mass and A 
a typical hadronic mass scale, the mass dependence of the QCD partition 
function is given by the effective partition function 

Z cS (M,V)= I dU exp(Ret/Str MUe l9/N f). (6) 

JUeG/H 

Here, M is the mass matrix, 9 is the vacuum angle and £ = | (qq)\. The integra- 
tion is over the Goldstone manifold G/H which is determined by the scheme 
of chiral symmetry breaking. In QCD with three colors and Nt fundamental 
fermions, G/H = SU(N f ). 

However, QCD is not the only theory that can be reduced to the effective 
partition function (^). In particular, chiral random matrix theories can be 
reduced to this partition functional. This leads us to the notion of universality: 
a property will be called universal if it is the same for theories with the same 
low energy effective theory. This allows us to calculate such quantities for the 
simplest possible theory, for example, for a chiral random matrix theory. In 
particular, properties that are determined by the effective partition function 
are trivially universal. As an example, I mention the Leutwyler-Smilga sum 
rules. Let us for the moment consider the sector of zero topological charge. 
Then, by expanding both sides of the equality (valid in the range @) 

Z?™(m,V) _ Zf =a (m,V) 



Z^(m^0,V) Zl%{m = Q,V) 

in powers of m, and performing the group integrals we obtain an infinite family 
of sum rules. The simplest sum rule is given byE3 

(8) 
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The generalization of this formula to a sector with topological charge v and 
different coset spaces {G/H = SU(2N f )/Sp(2N f ) for SU(2) color with fun- 
damental fermions, G/H = SU(Nf)/ SO(Nf) for adjoint fermions with gauge 
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group SU(N C ), N c > 2) is given by 

^ £, ( )* = 4 (M + (dim(cosct) + 1)/^) ' (9) 
The sum in @ can be replaced by an integral over the spectral density 

^E(w) =s 2 f §^L>. (10) 



V" 2 ^ \A?/ y u 2 VZ ry VY,' 

A fc >0 N fc ' v J ° 

where we have introduced the microscopic variable u = WE. The combination 

will be called the microscopic spectral density. It is constrained by the effective 
partition function but not determined by it. Our conjecture is that, in the 
thermodynamic limit, it is universal function that can be calculated with the 
help of chiral random matrix theory. 

A more intuitive interpretation of ps(u) is obtained if we start from the 
Banks-Casher relation and the fact that chiral symmetry is broken. The eigen- 
values near zero virtuality are roughly equally spaced, with spacing given by 

AX =m=^- (12) 

and smallest nonzero eigenvalue of the Dirac operator is equal to 

Aniin ~ r—;. (13) 

The microscopic variable u parameterizes the spectral density on the scale of 
the average level spacing, and the microscopic spectral density is obtained by 
magnifying the region around A = by a factor V . 

An observable directly related to the microscopic spectral density is the 
valence quark mass dependence of the chiral condensate EJ which can be ex- 
pressed as 

S(m) f°° 2mVZ , . , 

It can be obtained from random matrix theory in the range This range 
(||) can also be written as 



<A, (15) 
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which might be called the mesoscopic window of QCD. 

To summarize this section, we conclude that the eigenvalues of the QCD 
Dirac operator are strongly correlated for the following reasons: 

• The Leutwyler-Smilga sum-rules. 

• The eigenvalue spacing near zero AA ~ 1/V instead of AA - l/V 1 / 4 for 
a noninteracting system. 

• If the eigenvalues where uncorrelated, then, in the chiral limit p(X) ~ 
X 2N f , and chiral symmetry would not be broken. 

3 Universal level correlations of quantum spectra 

Quantum spectra of complex systems have been investigated in great detail 
both experimentally and theoretically. Already, in the fifties it was observed 
that the spectral correlations of nuclear resonances E3 do not depend on the 
details of the system and can be obtained from random matrix theory. More 
recently, it has been shown that the essential ingredient to obtainjandom 
matrix spectra is that the corresponding classical system is chaotic 13. Even 
systems with only two degrees of freedom show spectral correlations that are 
in perfect agreement with the random matrix prediction. It was already noted 
early on that the average spectral density is not given by random matrix the- 
ory. Typically, for a physical system, it is a monotonously increasing function 
whereas for random matrix theory it is a semicircle. To deal with this problem, 
a first important hypothesis was introduced, namely, that we have a separa- 
tion of scales with regards to the average spectral density and the microscopic 
spectral fluctuations. The second hypothesis is that the microscopic fluctua- 
tions are given by the random matrix ensemble with same symmetries as the 
underlying theory. In particular, anti-unitary symmetries play an important 
role in identifying the correct RMT. 

In practice, the average spectral density is folded out by transforming the 
original spectrum to a new spectrum with average spectral density equal 
to one. Below we will use the following statistics to measure the spectral cor- 
relations of the unfolded level sequence: the nearest neighbor spacing distribu- 
tion, P{S), the number statistics and the A3-statistic. The number statistics 
are obtained by counting the number of eigenvalues rik in consecutive intervals 
of length n. The number variance is defined as T,2(n) = var{rik}, and the As- 
statistic is obtained from a convolution of £2(71) with a 'smqathening'-kernel. 

The above statistics have been obtained analytically (seeliU) for the invari- 
ant random matrix ensembles defined by a matrix of independently distributed 
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gaussian matrix elements, consistent with the hermiticity of the matrix. De- 
pending on the anti-unitary symmetry the matrix elements are real (Gaussian 
Orthogonal Ensemble (GOE), (3 = 1), complex (Gaussian Unitary Ensem- 
ble (GUE), (3 = 2) or quaternion real (Gaussian Symplectic Ensemble (GSE), 
(3 = 4). The most notable property is the stiffness of RMT spectra: the number 
variance behaves asymptotically as £2(71) ~ log(n) instead of Y 2 (n) ~ n 
for independently distributed eigenvalues. 



4 Chiral random matrix theory 

Chiral random matrix theories are theories with the global symmetries of the 
QCD Dirac operator, but otherwise independently Gaussian distributed ran- 
dom matrix elements. n Xhe chiral random matrix model that obeys these con- 



ditions is defined byt 

r Nf NY? 13 
Zl= I DW \\ det(X> + m f ) exp( -tL^W^W), (16) 



where 



/ iW + i(n T -arg(P))+ti\ , , 

\ iW^ + i(tt T - arg(P)) + /z) J ' y ' 

and W is a rectangular n x m matrix with jz = n — m and N = n + m. An 
early version oLthis model can be found inE3. The temperature dependence 
introduced in BEI is included in the form of Matsubara frequencies in the 
diagonal matrix £It- The chemical potential, /1, introduced inE3 results in 
a non-Hermitean Dirac operator with eigenvalues that are scattered in the 
complex plane. The term. arg(P), the argument of the Polyakov loop was 
introduced by StephanovEil in order to explain the dependence of the critical 
temperature on the expectation value of the Zn -phase c3. In this paper we 
only discuss models with fi = and temperature dependence with only the 
lowest Matsubara frequencies included. This partition function reproduces the 
following symmetries of QCD: 



The Ua(^) symmetry. All nonzero eigenvalues of the random matrix 
Dirac operator occur in pairs ±A. 

The topological structure of the QCD partition function. The matrix 
T> has exactly \u\ = \n — m\ zero eigenvalues. This identifies v as the 
topological sector of the model. 
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• The flavor symmetry is the same as in QCD. For j3 = 2 it is SU(Nf) x 
SU(Nf). For (3 = 1 it is SU(2N f ) and for /3 = 4 it is SU(N f ). 

• The chiral symmetry is broken spontaneously for two or more flavors 
according to the pattern SU(N f ) x SU(N f )/SU{N f ), SU(2N f )/ Sp(N f ) 

Sd SU(Nf)/0(Nf) for /3 = 2, 1 and 4, respectively the same as in QCD 
. The chiral condensate satisfies the Banks-Casher relation 

Z = d m AogZ= lim ^P-. (18) 

./V— >OG iV 

Therefore, the parameter S in the random matrix model is identified as 
the chiral condensate and N as the (dimensionless) volume of space time. 

• The anti-unitary symmetries. For three and more colors with fundamen- 
tal fcrmions the Dirac operator has no anti-unitary symmetries, and, in 
general, the matrix elements of the Dirac operator are complex. The 
matrix elements of W of the corresponding random matrix ensemble are 
chosen arbitrary complex as well ((3 = 2). For N c = 2 and fundamental 
fermions the Dirac operator satisfies 

[Ct 2 K, i 1 D] = 0, (19) 

where C is the charge conjugation matrix and K is the complex con- 
jugation operator. Because, (CT2K) 2 = 1, the matrix elements of the 
Dirac operator can always be chosen real, and the corresponding random 
matrix ensemble is defined with real matrix elements ((3 = 1). For two 
or more colors with fermions in the adjoint representation ijD has the 
symmetry 

[CK, ijD] = 0, (20) 

but now (CK) 2 = — 1, which allows us to rearrange the matrix elements 
of the Dirac operator into real quaternions. The matrix elements of W 
of the corresponding random matrix ensemble are chosen quaternion real 
as well (13 = A). 

5 Universal spectral correlations 

As already stated before, we have to distinguish spectral correlations near 
zero and spectral correlations in the bulk of the spectrum. The latter have the 
advantage that, under the assumption of spectral ergodicity, they can be calcu- 
lated by a spectral average rather than an ensemble average. This requires no 
more than a few statistically independent gauge field configurations, whereas 
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the calculation of the microscopic spectral density requires a prohibitively large 
number of independent gauge field configurations. 

We will present three pieces of evidence that spectra of the lattice QCD 
Dirac operator show universal correlations on the finest scale. First, we will 
show that spectral correlations in the bulk of the spectrum are given by the 
invariant random matrix ensembles. Second, we will show that the valence 
quark mass dependence in the mesoscopic range of QCD can be obtained from 
the microscopic spectral density. Third, we will present random matrix results 
showing that that microscopic spectral density is insensitive to the probability 
distribution of the matrix elements and to the temperature as introduced in 
©• 

5. 1 Correlations in the bulk of lattice spectra 

Recently, conrolete Dirac spectra for reasonably large lattices were obtained 
by Kalkreuter E3 both for Kogut-Susskind fermions (with four flavors) and 
Wilson fermions (with two flavors) with SU (2)-color. It is important that all 
eigenvalues were obtained. Because the sum of the squares of the eigenvalues 
satisfy rigorous sum-rules, this allows a very stringent test on the accuracy of 
the calculation. 

In the case of the SU (2) color group, the anti- unitary symmetries of Kogut- 
Susskind fermions are different from the continuum theory. We have have 

[t 2 K,D ks }=0 with (t 2 K) 2 = -1, (21) 

(K is the complex conjugation operator) whereas for Wilson fermions in the 
form considered by Kalkreuter we have 

[7 5 CXr 2 ,757^ Wilson ] = with( 75 CXr 2 ) 2 = 1. (22) 

According to random matrix theory0, if the square of the anti- unitary operator 
is —1 the spectral correlations are given by the GSE, whereas, if the square is 
1, the correlations are given by the GOE. 

We have analyzed the spectra of a score of lattice QCD configurations. 
In Fig. 1 we show the results for the correlations of the Dirac eigenvalues for 
Kogut Susskind fermions (right) on a 12 4 lattice and Wilson fermions (left) on 
a 8 3 x 12 lattice for parameters as indicated in labels of the figure. We show 
the number variance, ^(n), the A3 statistic and the nearest neighbor spac- 
ing distribution P(S). The points represent the lattice-data and the random 
matrix results are given by the curves. The corresponding RMT is given in 
the label of the figure. We find perfect agreement with the above predictions. 
In addition, we have also confirmed that the correlations are stationary over 
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Figure 1: Spectral correlations of Dirac eigenvalues for Wilson fermions (left) and Kogut- 
Susskind fermions (right). The number variance, E2(n), the A3-statistic and the nearest 
neighbor spacing distribution are shown in the upper, middle and lower figure, respectively. 
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Figure 2: The valence quark mass dependence of the chiral condensate E v (m) plotted as 
T, v (m) /£ versus mNH. The dots and squares represent lattice results by the Columbia 
group [11] for values of /3 as indicated in the label of the figure. 

the spectrum. Even if we consider only the first 100 or 200 eigenvalues, no 
deviations from RMT correlations are found. 

5.2 The valence quark mass dependence of the chiral condensate. 

Using standard methods of random matrix theory it is possible to obtain an- 
alytically the microscopic spectral density. The result for Nf flavors in the 
sector of topological charge v isu 

Ps(u) = | {Jnj+u{u) - JN f +u+l(u)J Nf +„-l{u)j . (23) 

Using ( jl4| ) _ssje obtain the valence quark mass dependence by an elementary 
integrations 

— — = x{I Nf+u {x)K Nf+ v{x) + I Nf+v+ i{x)K Nf+v -i{x)). (24) 

Asymptotically, for large x = mUS, ~ S(l — (iVj + is)/x) and for x — > 0, 

E(x) - xY,/2(N f + v) for N f + v ^ and T,(x) xElog(x) for N f ± 

v = 0. The valence quark mass was calculated by the Columbia group liil 
for unquenched Kogut-Susskind fermions on a 16 3 x 4 lattice with m sca a = 
0.01, two flavors and three colors. In Fig. 2 we plot their results in terms of 
the dimensionless quantities suggested by the valence quark mass dependence 
(|J), namely as E(x)/E versus x. We observe that the lattice data fall on a 
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universal curve given by (^4|) for Nf + v = 0. The reason that Nf — is that 
Aminfl ~ l.Oe — 4 -C m sca a, so that the calculation is quenched on the scale 
of the smallest eigenvalues. It is also not surprising that v — because the 
fermionic zero modes are completely mixed with the rest of the states. For 
larger x the lattice data deviate from the universal curve. This was to be 
expected because in terms of x the domain of validity is given by 

x « X^-. (25) 
5.3 Universality in random matrix theory 

The microscopic spectral density was calculated for the simplest possible ran- 
dom matrix ensemble. Recently, the microscopic spectral density for Nf — 
v =3 was calculated for ensembles with a nongaussian probability distribu- 
tkmH. The authors found that the addition of a polynomial in tr(W^W) to 
the exponent of the probability distribution had no effect on the microscopic 
spectral density. We studied the random matrix model ( |l7| ) with only the 
lowest Matsubara frequencies includedlld. It can be shown that it is sufficient 
to include only the positive frequency. This model has a critical temperature, 
T c , above which p(0) = 0. In spite of the fact that the average spectral density 
changes drasticall^Jaetween T — and T — T c (it is given by the solution of 
a cubic equation we find that the microscopic spectral density remains 

unaffected. 

Numerically, we studied a variety other ensembles. In particular, I men- 
tion the chiral Cauchy ensemble, with a with a diverging second moment of 
the matrix elements. Nevertheless, we find the same universal microscopic 
spectral density. Finally, we wish to mention that according to general uni- 
versality arguments in random matrix theory, the correlations in the bulk of 
the spectrum of the chiral random matrix ensembles are given by the invariant 
random matrix ensembles a. Therefore, we find that the microscopic spectral 
density and the correlations in the bulk of the spectrum are mutually inclusive 
and both arc strongly universal. 

6 Conclusions 

We have argued that the spectral fluctuations on the finest scale are deter- 
mined by the global symmetries of the QCD partition function, and can be 
obtained from chiral random matrix theories. In particular, the microscopic 
spectral density determines the valence quark mass dependence of the con- 
densate in the mesoscopic range of QCD. In this range we find perfect agree- 
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ment with recent lattice QCD data. We have also observed that correlations 
in the bulk of the spectrum are given by the invariant random matrix en- 
sembles. This implies that the spectral fluctuations are strongly suppressed 
E 2 (n) ~ log(n)/-7r 2 , instead of n for uncorrelated eigenvalues. To some extent, 
QCD is self-quenching. 
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